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Beta-decay rates for spherical neutron-rich r-process waiting-point nuclei are calculated within a 
fully self-consistent Quasiparticle Random-Phase Approximation, formulated in the The Hartree- 
Fock-Bogolyubov canonical single-particle basis. The same Skyrme force is used everywhere in 
the calculation except in the proton-neutron particle-particle channel, where a finite-range force is 
consistently employed. In all but the heaviest nuclei, the resulting half-lives are usually shorter 
by factors of 2 to 5 than those of calculations that ignore the proton-neutron particle-particle 
interaction. The shorter half-lives alter predictions for the abundance distribution of r-process 
elements and for the time it takes to synthesize them. 

PACS numbers; 21.30.Fe, 21.60. Jz, 26.30.-fk, 23.40.-s, 23.40.Hc 



I. INTRODUCTION 

The astrophysical r-process [0-^ , which creates about 
half of all nuclei with ^ > 70, proceeds through very 
neutron-rich and unstable isotopes produced by stellar 
explosions or other violent events. The ultimate abun- 
dance of any stable nuclide depends strongly on the beta- 
decay half lives of its neutron-rich progenitors. The so- 
lar elemental abundance distribution shows peaks near 
y4.=80, 130, and 195, corresponding to progenitors with 
closed neutron shells (iV=50, 82, and 126). These rel- 
atively long-lived nuclei not only define the abundance 
peaks but also restrict the amount of heavier material 
that is synthesized. Understanding the important fea- 
tures of the r-process therefore requires knowledge of life- 
times of closed-shell semimagic nuclei far from stability. 
Of course beta decay is only one of the processes that 
contribute to r-process abundances; neutron capture and 
photodisintegration also play important roles, as do the 
temperature, density and initial neutron to seed ratio in 
the explosive environment. But these other aspects are 
beyond the scope of this article, which focuses on the 
crucial question of how to calculate beta decay far from 
stability. At present, the very neutron-rich nuclei are 
out of experimental reach, and theory provides the only 
handle on their decay rates. 

Precise theoretical estimates of beta-decay rates are 
hard to make. Most of the strength associated with the 
/3~-decay operator arj^ lies in the Gamow- Teller (GT) 
resonance, well above decay threshold. The strength that 
actually contributes to beta decay is the small low-energy 
tail of the GT distribution. Calculated beta-decay rates 
can therefore vary over a wide range without coming into 
conflict with sum rules, which for other processes help 
reduce theoretical uncertainty. In addition, beta-decay 



lifetimes depend sensitively on nuclear binding energies, 
and small errors in calculated Q^-values can cause large 
errors in predicted decay rates. These problems are com- 
plicated enough to demand a coherent and systematic 
approach to beta decay. Here we make a first attempt 
at a completely self-consistent calculation. Our goal, in 
spite of the difficulties, is a reliable estimate of beta-decay 
rates in important even-even semimagic nuclei lying on 
the r-process path. 

Special tools are needed to describe transitions to low- 
lying excited states in weakly bound nuclei. Although 
large-scale shell model calculations successfully describe 
the GT strength distribution in medium-mass nuclei close 
to the valley of beta stability [Rji], large configuration 
spaces and difficulties with the continuum |S] make the 
approach hard to apply along the r-process path. The 
continuum random phase approximation ||lO,ni| is often 
useful, but in very-weakly bound nuclei pairing is im- 
portant and a quasiparticle random phase approxima- 
tion (QRPA) based on coordinate-space Hartree-Fock- 
Bogolyubov (HFB) theory is required. Surprisingly lit- 
tle work has been done along these lines. Much more 
common are global (in that they attempt to describe all 
nuclei in the same framework) but non self-consistent cal- 
culations li~ 



that substitute the Strutinsky method 
BCS for HFB, approximate the continuum by bound 
or quasibound orbits, and use a schematic interaction 
K(cricr2)(TiT2) in the QRPA. The work of Ref. Q, which 
successfully reproduced the half-lives of the nickel iso- 
topes and of ^^^Sn in the Tamm-Dancoff approximation, 
used self-consistent single-particle energies and orbits but 
retained the traditional schematic residual interaction. 
But like the global calculations, this work did not in- 
clude a proton-neutron {pn) particle-particle interaction 
in the QRPA. 

Borzov et al. ||l^,|l8| did use a more self-consistent 



method in restricted regions of the nuchde chart. Their 
starting point was an energy-density functional opti- 
mized for the regions they considered, and a reahstic 
interaction (inchiding a zero-range pn particle-particle 
component) in the QRPA. Their energy functional was 
spin- independent, however, and so they did not obtain 
the residual interaction from the second derivative of 
the functional, as required if the QRPA is to repre- 
sent the small amplitude limit of time-dependent HFB 
(this is what we mean by self-consistency in an HFB + 
QRPA calculation). Instead they chose a phenomeno- 
logical spin-dependent residual interaction p^ . In addi- 
tion, they neglected positive-energy single-particle states 
in the HFB calculation, though they did include the 
entire particle-hole continuum in the QRPA. In any 
event, Borzov et al. |T9,Bfl] have now abandoned attempts 
at fully self-consistent calculations in favor of the Ex- 
tended Thomas-Fermi with Strutinsky Integral (ETFSI) 
approach. No previous papers, in sum, have used the 
same interaction in the mean-field calculation and the 
QRPA, and most have neglected the residual pn particle- 
particle interaction. Neither has any previous work prop- 
erly included the effects of the low-energy continuum on 
pairing. 

Here we attempt to do better: we use the same inter- 
action in the HFB and QRPA calculations so as to pre- 
serve the small- amplitude limit of time-dependent HFB. 
In the pn particle-particle channel of the QRPA, we use 
a finite-range interaction to soften the divergences in- 
duced by a delta function, but will show that this choice 
does not spoil self-consistency. By adjusting one param- 
eter, the strength of the pn particle-particle interaction 
(which we will often refer to as T=0 pairing), we repro- 
duce measured decay rates near the closed neutron shells, 
and then make predictions for rates further from stabil- 
ity. For now we restrict our study to semi-magic nuclei 
with closed neutron shells and their neighbors; besides 
being the most important for the r-process these nuclei 
are conveniently spherical. 

Our paper is organized as follows: In Sec. O we re- 
view the the general formalism for calculating beta de- 
cay, briefly present coordinate-space HFB theory and 
the QRPA, and then show how to combine them in a 
self-consistent way. Section IH discusses the problem of 



choosing a good effective interaction. Section IV contains 
the results of numerical calculations and an r-process sim- 
ulation, and Sec. ^ presents conclusions. 



II. FRAMEWORK 

A. Calculation of half-lives 

The rate for allowed Gamow- Teller decay of an even- 
even nucleus is given by 
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where the index n labels particles, lO"*") is the initial 
ground state and the |1+) are the final states (We use 
units in which c=l.) The sum over m runs over all 
final 1+ states with an excitation energy smaller than 
the Qp- value. We will set the weak axial nucleon cou- 
pling constant qa to 1 rather than its actual value of 1.26 
to account for the near universal quenching of isovector 
spin matrix elements [ pi[ in nuclei. The kinematic factor 
dnm/dE is the density of final (e~, De) states. It can be 
written as 
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where Ei is the ground state energy of the even-even 
initial nucleus with Z protons and N neutrons, E^+ is 
the energy of the mth excited l"*" state in the final nucleus 
with Z + 1 protons and N — 1 neutrons, and F{Z,Ee) 
corrects the phase-space factor for the nuclear charge and 
finite nuclear size, both of which affect the electron wave 
function. 

We would like an approximation for Ei ~ E^+ in Eq. 
(0) that does not require an explicit calculation of the 
Qp- value. To obtain one we express the Qp- value in 
terms of nuclear ground-state binding energies: 
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where AM„_//=0. 78227 MeV is the mass difference be- 
tween the neutron and the hydrogen atom. The ground- 
state binding energy of an odd-odd final nucleus, in the 
independent quasiparticle approximation, is given by 

Bg.s.{Z -f 1, A^ — 1) « i?g.s.(Z, N) + \p — Xn + E2qp^ovjost 

(4) 

where Xp=dE / dZ and Xn=dE / dN are the proton and 
neutron Fermi energies in the initial nucleus {Z, N), and 
E2qp,iowcst is the energy of the lowest two-quasiparticle 
excitation (with respect to the initial nucleus) corrected 
for the residual pn interaction between the valence quasi- 
particles in the odd-odd system [ P2| . The excitation en- 
ergies of the 1'^ states with respect to the final ground 
state are 
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where i^QRPA is the QRPA phonon energy (see Sec. ttlQ). 
It follows from Eqs. (||) and (^ that the energy released in 
the transition from the ground state of the initial nucleus 
to a l"*" state in the final nucleus is 



E,-E^+^=Qp_{Z.N)-El^ 

~ A„ - Ap + AMn-H - £^QRPA • (6) 

This expression allows us to avoid calculating ground- 
state masses of the final nuclei. 

B. Skyrme-Hartree-Fock-Bogoliubov Method 

Our first-order description of the even-even ground 
state and quasiparticle excitations is based on a self- 
consistent mean field. Several schemes are currently 
used in mean-field calculations for heavy nuclei, among 
them the non-relativistic HF and HFB methods with 
the Skyrme interaction [E3[ , the HFB approach with the 
Gogny force [g4|, and the relativistic mean-field model 
p5| , p6[ . All these provide a parameterized energy func- 
tional that describes nuclear properties throughout the 
chart of nuclides. But in neutron-rich nuclei near the 
drip-line, where the Fermi energy and pairing potential 
are close in magnitude, pairing cannot be treated as a 
small correction to the HF solution. Instead the physics 
underlying the HF and BCS approaches must be incor- 
porated into a single variational principle, through HFB 
theory. Our version of HFB is formulated in coordinate 
space and described in detail in Refs. [ p7p^ ]; here we dis- 
cuss only the key ingredients needed for calculating beta 
decay. 

Our effective interaction comes from the Skyrme 
energy-density functional £. The functional can be di- 
vided into two parts: 
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The particle-hole part £ph depends primarily on the 
density matrix pq{r(7,r'a'), while the pairing (particle- 
particle) part fpair depends primarily on the pairing den- 
sity matrix Xq(rcr, rV) (the index q denotes protons or 
neutrons). The coupling between p and x comes from 
the density-dependent terms of the Skyrme and/or pair- 
ing interaction. Because the nuclei we consider are so 
neutron-rich, we shall neglect pn pairing in the HFB 
ground state (see Sec. EH); as a result our HFB wave 
function will be a product of proton and neutron wave 
functions, and matrices of p and x "^iH be block diago- 
nal in q. Expressions defining the Skyrme and pairing 
energy functionals have been presented many times (see 
e.g. Ref s. [P9|-p2[ ) and so are only briefiy discussed here 
in Sees. |III A| and |lll4 

The particle-hole, hq{ra,r'a'), and particle-particle, 
Aq{ra,r'a'), mean fields are defined as the first deriva- 
tives of the energy functional with respect to the corre- 
sponding densities: 
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The residual interactio ns th at enter the QRPA equa- 
tions (discussed in Sec. ilC) are the corresponding sec- 
ond derivatives of the energy functional. Rearrangement 
terms, which result from the density dependence of the 
energy functional are therefore included in the both the 
HFB and the QRPA. 

Variation of the energy functional (R) with respect to 
the densities leads to the coordinate-space HFB equa- 
tions for protons and neutrons p7|, 
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where a enumerates the HFB quasiparticle eigen- 
states. The two-component quasiparticle wave func- 
tions 0i(a,r'(T') and (/)2(a, r'cr') self-consistently define 
the densities: 



Pq(rcr,r'cr') = ^(/)2(a,r,cr) (t)*{a,r' ,a') 
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The Lagrange multipliers A^ — the Fermi energies of the 
neutrons and protons — are fixed by the particle-number 
conditions 
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By working directly in the coordinate space we are able 
to properly include unbound states, which, as we have re- 
marked, become important near the neutron drip line. A 
particular virtue of our approach is the accurate repre- 
sentation of the canonical single-particle basis, consisting 
of eigenstates V'g/j i'f''^) oi the density matrix: 
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Pq{ra, r'a')^q^{r'a') = v^ iPq^{ra). (14) 



Canonical states form an infinite, discrete, and complete 
set of localized wave functions ||2^,|2^ ; they describe both 
the bound states and the positive-energy single-particle 
continuum. The canonical basis is well defined for all 
particle-bound nuclei, no matter how close they are to 
the drip line. We use it in the next section to formulate 
the QRPA. 



C. Proton- neutron QRPA 

Although one can derive coordinate-space QRPA equa- 
tions for the generalized density matrix [ pO| , pT|p^ ], here 
we stick with the older representation in terms of a dis- 
crete single-particle basis. The canonical HFB basis is 



ideal for this purpose because it simplifies the HFB equa- 
tions, so that HFB + QRPA can be formulated in com- 
plete parallel with ECS + QRPA (described, e.g., in Ref. 
|p3|). The pn QRPA equations take the form: 
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with matrices A and B defined as 
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Here p, p' , and n, n' denote proton and neutron quasipar- 
ticle canonical states, V is the pn particle-hole interaction 
in the 1+ channel, obtained from the second derivative 
of energy functional £ with respect to the proton and 
neutron densities, and V is the corresponding particle- 
particle interaction, obtained from the second derivative 
with respect to the pairing densities. The X's and y's are 
the amplitudes for exciting two-quasiparticle and two- 
quasihole states from the correlated vacuum. The oc- 
cupation amplitudes Vi are the eigenvalues of the den- 
sity matrix (p^). Since the canonical HFB basis does 
not diagonalize the HFB Hamiltonian ([l(i|), the one- 
quasiparticle terms in the Eq. (|lj) have off-diagonal ma- 
trix elements Eij . The presence of these terms is the only 
formal difference between our QRPA equations and those 
based on the BCS approximation. The localized canoni- 
cal wave functions, however, are more realistic than the 
single-quasiparticle states used to mock up the contin- 
uum in the BCS approximation. 

To actually solve the matrix QRPA equations (^ we 
have to truncate the canonical basis. The diagonal ma- 
trix elements of the one-quasiparticle Hamiltonian pro- 
vide a convenient measure of excitation, and we include 
only those canonical states for which these matrix ele- 
ments are less than a cut-off value. (For the choice of 
the cut-off energy, see Sec. IH.) Occupation probabili- 
ties of canonical states usually decrease much faster with 
excitation energy than in the BCS approximation P8[ , 
allowing us to work with matrices of manageable size. 



III. INTERACTIONS 

Different channels of the effective interaction have dif- 
ferent effects on the beta-strength function. Roughly 
speaking, the like-particle and pn particle-hole interac- 
tions determine the single-particle spectrum. The like- 
particle pairing interaction smears the Fermi surface 
and changes the elementary excitations from particles to 



quasiparticles. Without further many-body correlations, 
the GT resonance is not collective and too much of the 
beta strength is located at low energies. The pn particle- 
hole force, treated in the QRPA, solves the problem by 
pushing GT strength up and into the resonance, leav- 
ing the low-lying spectrum depleted. This depletion is so 
great, however, that the resulting lifetimes can become 
too long; the half-lives of Refs iQ, for example, usually 
exceed measured half-lives. The pn particle-particle force 
cures this problem by pulling some strength back down. 
Interestingly, this part of the force is also necessary to 
describe /S"*" decay in proton-rich nuclei, but there it de- 
creases the low-lying strength. 

To be consistent, one should use a single interaction 
in all channels and in every step of calculation. This 
means, for example, that the same pn particle-hole en- 
ergy functional that determines single-particle properties 
in the HFB calculation must also be used in the QRPA. 
The same is true for the pairing (particle-particle) forces, 
both in the like-particle and pn channels. The constraint 
is not as tight as in the particle-hole channel, however, 
because the proton and neutron Fermi surfaces are so 
far apart in neutron-rich nuclei that pn pairing corre- 
lations are neglibigle in the HFB. Furthermore, to the 
extent that the T—l pairing force can be approximated 
by a delta function, it does not affect the J'^=l+ states 
obtained in the QRPA. Even with a finite range force, 
the effects are negligible. The pn component of the T=l 
pairing/particle-particle interaction can therefore be ne- 
glected everywhere. In addition, one is free to choose the 
T=0 pairing component solely on the basis of its effects 
in the QRPA, since it has no like-particle component and 
does nothing at the mean field level unless Nk,Z. 

For these reasons we organize the discussion in this sec- 
tion as follows: First we describe the particle-hole inter- 
action in detail; it largely determines both single-particle 
properties and the collectivity of the GT resonance. We 
then briefiy discuss the like-particle pairing force, which 
plays a role only in the HFB calculation. Finally, we 
describe the pn (T=0) pairing interaction. This last in- 
gredient appears only in the QRPA but, as mentioned al- 
ready, is crucial for reducing calculated lifetimes to values 
that are consistent with experiment. 



A. The particle-hole interaction 

The pn particle-hole interaction is responsible for the 
main features of the GT distribution. For our calcu- 
lations to make sense, we need to find an interaction 
that reproduces the distribution reasonably, at both high 
and low energies. Before doing this, we have to discuss 
the general form of the particle-hole interaction in the 
Skyrme framework. 

For even-even nuclei, the particle-hole Skyrme energy 
functional of Eq. (0) can be written as 
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i.e. as the sum of a kinetic-energy functional £'kin, an ef- 
fective strong-interaction functional fsk, and a Coulomb 
functional £c- The functional are the spatial integrals 
of the corresponding local energy densities 7i, 



of the local density approximation [p7 38 1 , but has been 
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The kinetic-energy and Coulomb energy densities are 
given by 
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The Skyrme energy density can be split into pieces 7i^ 
and Wgk'^ that are bilinear in time-even and time-odd 
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densities, respectively [|0|3J]. Only Hg];™ affects the 
ground-state properties of even-even nuclei; it can be 
written as 
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The energy densities depend on the local matter density 
Pq, the kinetic density r^, and the spin-orbit current J^. 
Densities without a q index are total (isoscalar) densities, 
e.g. Pq = Pp+Pn- The parameters bi, b[, ci, and c'l are the 
time-even coupling constants of the Skyrme functional. 
The parameters all play different roles; the terms with 
63 and ^3, for example, determine the density-dependent 
parts of the interaction, and the 64, 64 and ci, c'l terms 
define the spin-orbit interaction. All the parameters are 
fit to a few key nuclear-structure data, e.g. total binding 
energies, charge radii, and surface thicknesses of nuclei 
in the valley of stability, and they generally reproduce 
ground-state properties of nuclei between ^^O and the 
heaviest elements [ p5| , p6[ . 

The energy density TigjJ^, the detailed form of which 
can be found in Ref. [g^, involves time-odd quantities: 
the momentum density, spin density, and vector part of 
kinetic-energy density. The time-odd terms contribute 
to the energy of polarized states, i.e. those with nonzero 
angular momentum, including the 1"*" states populated 
by beta decay. As a consequence, the distribution of GT 
strength depends primarily on this time-odd part of the 
energy density. Although most of the coefficients in TigjJ'^ 
are in principle independent parameters, they are fixed 
by the values of the parameters in Tig™" in the usual 
Skyrme ansatz |34]. The restriction is not in the spirit 



made implicitly in almost all prior work. The predictions 
of most existing Skyrme forces for excited states in odd- 
odd nuclei are therefore completely determined by fits 
of Hg™" to properties of even-even ground states. Not 
surprisingly, we find the GT spectra predicted by most 
of these forces to be well off the mark. 

Clearly, the proper approach to beta decay in the 
Skyrme framework is to fit coefficients in 'H'^ to prop- 
erties of nuclear states with non-zero angular momen- 
tum. That task is not so easy, however. Excited states 
take more computer time to explore than ground states 
because they require methods beyond mean-field theory 
(such as the QRPA). Furthermore, 'H'^^ contains many 
parameters, and fixing them all would be a major under- 
taking. For these reasons, we have chosen here to use 
the relations between time-even and time-odd densities 
imposed by the the traditional Skyrme- force ansatz, in 
spite of the drawbacks. (The ansatz, incidentally, relates 
the 6i, &^, Cl, and c'l of Eq. ( p2| ) to the usual Skyrme- force 
parameters ti and Xi, as spelled out in Ref. |3|].) We will 
explore the effects of independent time-odd densities in 
a future publication. 

Even with constraints on TigjJ'^, it is difficult to include 
excited-state properties in the data set used to fit Hg™". 
Instead we follow the approach advocated in Ref. I^^l and 
compare the effective strength of our spin- and isospin- 
flip interaction in nuclear matter with phenomcnological 
values or predictions of realistic calculations. In nuclear 
matter the interaction strength in the GT channel is sen- 
sitive to a single combination of the Skyrme parameters, 
the Landau-Migdal parameter g^ [^ : 
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where we have used the normalization factor Nq — 
2kFTn* /{ir^h ), with the Fermi momentum defined as kp 
= (Stt^ /5/2)^/'^. Table || shows predictions of represen- 
tative Skyrme interactions - SGII [|9|; SkM* @; SkP 
13; SLy4, SLy5 Q; SkM ||; and SkO' ^ - for prop- 
erties of symmetric nuclear matter at saturation density. 
While all the forces give about the same saturation den- 
sity Pa, energy per nucleon E/A, and incompressibility 
Koo, small but noticeable differences appear in the ef- 
fective mass m* and asymmetry coefficient ii'ooQ The 
Landau parameter /o is a combination of po, m* and 
Koo, and /g is related to po and to*; it is therefore not 
surprising that most forces predict similar values for the 
time-even quantities /o and fL and that those values 
are close the empirical ones ESI. But the parameters 



^Not all the properties in the Table 1 are predictions of the 
parameter sets; the effective mass of SkP, the incompressibil- 
ity of SLy4 and SLy5, and the asymmetry coefficient of SLy4, 
SLy5, and SkO' have been used as constraints in the fits. 



TABLE I. Properties of symmetric nuclear matter at saturation density po predicted by representative Skyrme interactions. 
EjA is the energy per nucleon, m* /m the effective mass, Koo the incompressibility, and asy-m the asymmetry coefficient. 
The remaining quantities are the Landau- Migdal parameters fo,f}ugo, and q'q [BQI; q'q controls the effective interaction in the 
spin-isospin channel. The empirical values were taken from Ref. |43] 



Force 


po 

[fm-^] 


E/A 

[MeV] 


m* /m 


K^ 


Q.syni 

[MeV] 


/o 


/o 


go 


90 


SGII 


0.158 


-15.58 


0.786 


214 


26.8 


-0.233 


0.728 


0.622 


0.934 


SkM* 


0.160 


-15.75 


0.789 


216 


30.0 


-0.229 


0.926 


0.325 


0.937 


SkP 


0.162 


-15.93 


1.000 


201 


30.0 


-0.102 


1.417 


-0.229 


0.062 


SLy4 


0.160 


-15.97 


0.694 


230 


32.0 


-0.276 


0.813 


1.385 


0.901 


SLy5 


0.161 


-15.98 


0.698 


230 


32.0 


-0.276 


0.814 


1.137 


-0.152 


SLy6 


0.159 


-15.92 


0.690 


230 


32.0 


-0.280 


0.803 


1.408 


0.899 


Skf4 


0.160 


-15.92 


0.650 


248 


29.5 


-0.273 


0.559 


1.768 


0.881 


SkO' 


0.160 


-15.75 


0.896 


222 


32.0 


-0.097 


1.328 


-1.612 


0.792 


Empirical 












0±0.3 


~1.6 


~0.4 


~1.80 



(7o and ^g, which act in time-odd channels, reflect spin- 
dependent components that are not constrained by stan- 
dard fits. As a result their values scatter within a wide 
range: -1.6 <go< +I.S and -0.2 < g^, < +1. Further- 
more, two otherwise very similar forces can give quite 
different values for these parameters. The interactions 
SLy4 and SLy5, for example, are identical for one small 
detail, the treatment of the J^ term during the fit. But 
despite the close agreement for most nuclear-matter pa- 
rameters, the differences in both g^ and g^ between the 
two are significant. 

We remarked above that predictions of existing Skyrme 
forces for GT distributions usually are not good. This 
fact is reflected in Table | by the values of g'^ , all of which 
are much smaller than the empirical value (7o(expt) « 1.8 
p3[ . In fact, we are not aware of a single Skyrme inter- 
action that gives g'g close to the empirical value and at 
the same time does an acceptable job with global nuclear 
properties and single-particle shell structure. State-of- 
the-art Skyrme forces tend to yield a g'^ of about 0.9, a 
value too small by a factor two. An important related 
point: the value 5q=0.503 reported for SGII in Ref. ||3^] 
does not correspond to the consistent application of any 
Skyrme force. The authors left out the last term (in 
brackets) in Eq. (E2) when fitting the parameters of their 
force, but included its effects in their calculation of ^g 
(and in their RPA calculations). To be consistent with 
their mean-field interaction, they should have omitted 
the term proportional to kj^ in Eq. (ESI). Doing so gives 
5g=0.93, the value in Table |. We will elaborate on this 
remark in a future paper. For now, what's important is 
that even calculated properly the value of g'g associated 
with this interaction, which was designed explicitly for 
GT resonances, is far too small. 

Of course g^ in symmetric nuclear matter at saturation 
density does not directly measure the strength of spin- 
and isospin-flip interactions in a finite nucleus. All effects 
due to the nuclear surface, finite nuclear volume, and ex- 
cess neutrons disappear in nuclear matter. Moreover, the 
GT distribution depends strongly on the single-particle 



spectrum as well as the residual interaction. Nonetheless, 
the parameter g'g summarizes the gross features of the 
distribution. This is illustrated in Figs. and 0, which 
display the summed GT strength calculated for Zr and 
^^^Cd. Here we use three Skyrme forces with different 
values of Landau-Migdal parameter: SkO' (5q=0.792), 
SkP (^1^=0.062), and SLy5 (g(,=-0.152). In general both 
the energy of the GT resonance and the strength in it in- 
crease with 5q; they are smallest for SLy5 and largest for 
SkO'. Forces with small values of g'g fail to concentrate 
enough strength in the resonance, leaving too much at 
low energy. 

Although none of the interactions have large enough 
5q when used self-consistently, some are better than oth- 
ers. Having opted for now not to increase gg by using 
energy functionals in which Hg^*^" and TYgJ^^ decoupled, 
we decided to use SkO' [Q, which has one of the larger 
values in Table m. The SkO' energy functional is more 
general than the original Skyrme functionals in the spin 
orbit channel; the extra generality manifests itself as a 
difference in the values of &4 and 64 ^^. Without this 
extended interaction it seems to be impossible to obtain 
a consistent description of spin-orbit splittings and other 
global observables unless the last term in brackets in Eq. 
(p3) is neglected |M. The extension remedies the prob- 
lem quite nicely, and in spite of the low value of gg, SkO' 
reproduces GT spectra [HS] measured in charge- exchange 
reactions fairly well. The predicted resonances, though 
usually a little too low in energy, contain about the cor- 
rect fraction of the strengthQ 



^This observation is important because most published 
Skyrme-RPA calculations do not address the problem of the 
GT strength distribution. For instance, Ref. |4q] proposes 
that the difference between the centroid energy of the GT 
resonance and the energy of the isobaric analog state be in- 
cluded as a constraint on effective interactions. Getting the 
right energy difference seems to be too weak a criterion; the 
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FIG. 1. Summed GT strength up to 20 MeV as a function 
of excitation energy for the closed-shell nucleus ^"Zr, calcu- 
lated with the SkO', SLy5, and SkP Skyrme forces. We also 
plot the measured strength reported in Ref. |59| The calcu- 
lated strength, as is customary, is multiplied by (1/1.26)^; the 
quenching corresponds to setting gA to 1.0 in our calculations 
of beta decay. 
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FIG. 2. Same as Fig. hi, but for the open-shell nucleus 
^^*Cd. No experimental data exist here. 



VhP 



194.6 MeV fm^ 



Vi, 



-186.7 MeV fm^ 



B. Like- particle pairing interaction 



C. Residual pn particle-particle interaction 



The T=l pairing interaction between like particles af- 
fects only the HFB part of our calculation. To incor- 
porate it we use a simple pairing energy functional that 
corresponds to a delta force, 



f, 



pair 



q—p^n 



Vi,,xlir) 



(24) 



where Xg(r) ~ So- Xgli^j f''! i"i '^) is the local pairing ten- 
sor. This interaction has vanishing matrix elements in 
the 1+ channel and therefore contributes nothing to the 
residual interaction in the pn QRPA. We adjust the 
strengths Vi^q in the HFB calculation to reproduce ex- 
perimental pairing gaps as explained in Ref. |17[ | , though 
unlike the authors of that paper we use different val- 
ues for protons and neutrons, and allow them to depend 
slightly on mass. This refinement appears to be neces- 
sary for a precise description of the beta decay rates. For 
light nuclei with A^ « 50 we adopt the values 



Vi, 



488.1 MeV fm^ 



Vi, 



-213.8 MeV fm-^ 



while for heavier nuclei with N oi 82 we use 



force SGII, which passes the test in Ref. Ma 
(?o=0.93 that is well below the empirical one. 



has a value of 



Although the GT resonance is built almost entirely of 
particle-hole excitations, the low-lying strength respon- 
sible for beta decay involves pn particle-particle correla- 
tions and is sensitive to the T=0 pairing interaction. The 
reason these correlations are important at low energies is 
that the proton orbitals near the Fermi surface are nei- 
ther completely empty nor completely occupied. They 
therefore can accept the additional particle created from 
occupied neutron orbitals by beta decay at the same time 
as they interact with those neutron orbitals through the 
T—0 pairing force. A level that is completely full, by con- 
trast, can interact with the occupied neutron orbitals but 
will not participate in beta decay, while one that is com- 
pletely empty can accept additional protons from beta 
decay but will experience no particle-particle interaction 
with the occupied neutron levels. 

Because T=0 pairing has no effect in our HFB calcula- 
tions, we can treat its strength as a free parameter in the 
QRPA. Our procedure is to fit that strength to known 
beta-decay lifetimes in regions of the nuclear chart near 
those that interest us. To keep matters simple, we re- 
strict ourselves to a density-independent force acting only 
in the S=l channel; QRPA calculations of double-beta 
decay [|3| and single-beta decay |^,|3[ have shown this 
component of the interaction to have the largest effect 
on low-lying strength. Unfortunately, although a simple 
delta force successfully describes like-particle pairing, it 
turns out to be inadequate here; the calculated half-lives 
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FIG. 3. Dependence of the half-life of ^^^Cd on the size 
of the canonical basis used in the QRPA, plotted versus the 
strength of the T=0 pairing interaction. The lower and upper 
panels show results obtained with the finite-range force of 
Eq. (Ea) and a delta force. The configuration space in QRPA 
is defined by the cut-off energy, Ccutoff as described in the text. 
The horizontal line shows the experimental half-life of ^^*Cd. 



diverge steadily as canonical single-quasiparticle states 
are added to the basis used in the QRPA. Actually, any 
purely attractive interaction suffers from the same prob- 
lem. The situation improves considerably, however, when 
we use (as a crude mock up of a microscopic G-matrix) a 
short-range repulsive Gaussian combined with a weaker 
longer-range attractive Gaussian: 



Fi2 = -Vc 



2 



9j 



e ''12/^j 115^1 x=o 



(25) 



where IIs=i,t=o projects onto states with S — 1 and 
T = 0. We take the ranges /ii=:1.2fm and /X2=0.7fm of 
the two Gaussians from the Gogny interaction |Q, and 
choose the relative strengths 51 = 1 and 32 =-2 so that the 
force is repulsive at small distances. The only remain- 
ing free parameter is Vq; the overall strength. Figure 
depicts the rate of convergence of the predicted half-life 
of "'^^^Cd as a function of Vq for the interaction ( p5[ ) and 
for a 6 interaction. The latter gives results that do not 
converge when we increase ecutoff, the upper limit on 
the (diagonal) energy of canonical states we include in 
the QRPA. The finite-range force behaves much better. 
The calculations reported below were carried out with 
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FIG. 4. Ratios of calculated-to-measured half-lives for 
four nuclei near A'' = 50 as a function of the T=0 paring 
strength. The solid horizontal line corresponds to equal val- 
ues for measured and calculated half-lives, the short-dashed 
horizontal line to a factor of 2 difference, and the dotted line 
to a factor of 5 difference. 
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FIG. 5. Same as figure 0, but for nuclei near N = 82. 



the two-Gaussian^ force in Eq. (|25| ) and with ecutoff=25 
MeV for states above the Fermi surfaces (we include all 
states below the Fermi surfaces). 



■^Other choices are possible. For instance, we have checked 
that the convergence can also be achieved with a simpli- 
fied particle-particle interaction of Skyrme type with io = 
6.16971^0 MeV fm^ fj = -2.5849^0 MeV fm^ xo = a;i = 
0.25, and all others equal to zero. 
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FIG. 6. Comparison of calculated half-lives for Zn, Cd, and Ni isotopes with (HFB+QRPA+SKO') and without 
(HFB+QRPA+SkO', Vb = 0) the residual particle-particle interaction, and results from Refs. @ (FRDM+QRPA), jl|] 
(ETFSI+QRPA), and Q (HFB+TDA+SkP, only for the Ni isotopes), with experimental values taken from Ref. Q where 
possible. For the nickel isotopes recent results from [pIIm] are shown as well. When predicted half-lives are larger than 100 s, 
the FDRM collaboration reports only this lower bound, which is marked here with arrows pointing up. The ETSFI collaboration 
reports half-lives only when the predicted deformation (32 is less that 0.1. 



To fit the strength Vqj we use recently measured half- 
lives of neutron-rich nuclei in regions where the r-process 
path comes closest to the valley of stability (e.g., around 
^^Ni) [|0|-|2)). Fig. I displays the ratio of calculated-to- 
experimental lifetimes for three zinc isotopes and *^Ge, 
all near A^=50. The four lines intersect near a ratio of 1, 
showing that the experimental lifetimes can all be repro- 
duced with a single value of the T—Q pairing strength, 
Vb=230 MeV. We use this interaction strength to predict 
the lifetimes of nuclei (in this region) that are still further 
from stability. 

Fig. 1^ shows the corresponding results for 126,128, isoq^j^ 
all near 7V=82. The fit, while not as good here, is still rea- 
sonable and gives a best value for Vb of about 170 MeV. 
At that value we slightly underestimate the lifetime of 
"'^^''Cd (by a factor of 0.6). For ^^^Cd the discrepancy is 
greater — we underestimate its lifetime by a factor of 5 
— but that nucleus has such a small Q^- -value that a 
small error in the strength distribution can have a large 
effect on the rate. We therefore adopt the value Vb=170 
for the entire region. The reason the value of Vb is so 
different here might be connected with the quality of the 
single-particle spectra and Q-values predicted by Skyrme 
forces, and with the simplicity of our T=Q pairing inter- 
action. We may be compensating for such deficiencies 
by changing Vb- The fact that we have to change the 
T=0 pairing strength by such a large amount when go- 
ing from A^=50 to N=%2 demonstrates the sensitivity of 
calculated beta decay rates to other parts of the effective 
interaction. 



In the N—\2Q region there are not enough experimen- 
tal data to support a fit of the T—Q pairing strength, so 
we use the same value (Vb=170MeV) as in the iV=82 
region. This means that our predictions around A'^=126 
are less reliable than in other regions (see, however, dis- 



cussion in Sec. |V|) 



IV. RESULTS 



A. Measured half-lives 



Figure ^ displays the results of our calculations for nu- 
clei near the neutron closed shells, alongside the QRPA 
results without T=0 pairing (Vb=0), the results of Refs. 
||l9| , |l6 53 L and measured half- lives iQ. (The predictions 
of Ref. |19{ appear only when the nucleus is thought to 
have small deformation. Our numbers in these nuclei 
are more suspect than elsewhere, since we ignore defor- 
mation in our calculation.) Because we separately adjust 
the T=0 pairing strength in two regions, it is not surpris- 
ing that our results are usually closer to experiment than 
those of global calculations in which parameter values 
are kept fixed. But the errors in the global calculations 
are systematic; they almost always overestimate the life- 
times. In the case of Ref. [0,p3|, at least, we attribute 
the problem to the neglect of the T—Q pairing. As the 
figure shows, the results of Ref. ||l2|,^ are much closer 
to ours in most of the nuclei when we turn that force off. 
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FIG. 7. Calculated half-lives of the neutron-rich Ni iso- 
topes as a function of the strength of the T—Q pairing inter- 
action. 



In the Ni isotopes our lifetimes are also too large by 
factors of 3-5. Part of the reason is the weak sensitivity 
of these lifetimes to Vq, which in turn is due to the Z—28 
and 7V=40 shell closures. The proton orbitals important 
for beta-decay are completely empty and therefore couple 
through T—0 pairing only to neutron orbitals that are at 
least partly empty. Because of the closed shell at iV=40, 
the lowest such neutron orbit is 59/2- But the 179/2 neu- 
trons can only decay into the proton g orbitals, which are 
far above the proton Fermi surface. Thus no low-energy 
strength can be moved by the particle-particle interaction 
and the curves in Fig. 0, which shows the calculated half- 
lives of the nickel isotopes versus Vq, are all flat. They 
do not turn down until after the value of Vq (230 MeV) 
that is suitable for the other nuclei in the region, and 
close to the point at which the QRPA collapses. Such 
behavior clearly means that SkO' is not optimal in this 
region; the SkP interaction used in the calculations of 
Rcf. |16| apparently better predicts single-particle prop- 
erties and Q-values. Those calculations use no particle- 
particle force but, as discussed above, none is required in 
the Ni isotopes. 

Our HFB -I- QRPA theory violates particle number 
conservation, with the result that pairing correlations ar- 
tificially break down at closed shells. The predictions for 
nickel would probably be better in a number-conserving 
version of our approach (for the general formalism, see 
Rcf. H). 



B. Closed-neutron-shell r-process w^aiting points 

The effects of T=0 pairing vary just as much in closed- 
neutron-shcU nuclei along the r-process path as they do 
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FIG. 8. The same as figure m, but for A'^ = 50 isotones. 



in the measured nuclei just discussed. In doubly magic 
nuclei like '^^Ni, ^^^Sn, and ^^^Zr, the T—0 pairing force is 
ineffective. On the other hand, when one takes away two 
protons from these nuclei, creating two holes in a high- 
spin orbit, the force has a large effect; the corresponding 
neutron orbit and its spin-orbit partner are not too far 
below the neutron Fermi surface and contain many neu- 
trons, which both interact with the protons at their Fermi 
surface and decay to fill the two proton holes. The effect 
of adding two protons to a closed shell is a bit smaller. 
The several orbits above the closed proton shell have 
lower spin or are far from the Fermi surface, and their 
contributions tend to cancel. These points are illustrated 
in Fig. H, which shows the dependence of the calculated 
half-lives of several iV = 50 isotones on Vq. The half-life 
of the doubly magic nucleus ^^Ni, of course, varies al- 
most not at all, and we probably overpredict its lifetime 
slightly just as in the other Ni isotopes. 

Figure ^ shows our predictions, together with those 
of other authors, for the half-lives of all the crucial even- 
even closed-neutron-shell nuclei along the r-process path. 
Our results agree fairly well with those of Ref. [n2| for the 
very proton-poor nuclei (with A^ = 50 and N = 82) but 
less well for larger Z. The trend is due to the closed pro- 
ton shells a,t Z = 20, 28, 40, and 50, where the particle- 
particle force has little effect. Between these magic num- 
bers, and particularly just below them (e.g., in ^^Fe), the 
differences can be large. To demonstrate again that they 
are due to T—0 pairing, we plot results once more with 
that component of the force switched off (Vo=0), a step 
that brings our results into agreement with those of Ref. 
||l2| in nearly all nuclei with iV =50 or 82. 

As discussed in Sect. Ill C , there are no experimen- 
tal data with which to fix Vq near iV=126. The lack of 
closed shells in this region suggests that half-lives will de- 
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pend strongly on Vq. Our results with and without T=0 
pairing, however, show that this is not the case. Even 
if we used much a smaller value of Vq, by extrapolating 
the drop in that parameter between A^ = 50 and 82, the 
lifetimes would not change appreciably. In these heavy 
systems our results agree well with those of Ref. ||l^ . 



C. Consequences for nucleosynthesis 

The closed-neutron-shell nuclei are instrumental in set- 
ting abundances produced in the r-process; new predic- 
tions for their half-lives will have an effect on the results 
of r-process simulations. For iV=50 and 82 our half- 
lives are usually shorter than the commonly employed 
half-lives of Ref. |l2|, and longer for N—12Q Replacing 
those lifetimes with ours should therefore produce smaller 
A « 80 and 130 abundance peaks, and a larger A « 195 
peak. 

Without extending our calculations to other nuclei in 
the r-process network, however, we cannot draw quan- 
titative conclusions from a simulation. Accordingly, we 
carry out only one simple r-process simulation here, com- 
paring final abundance distributions obtained from the 
beta decay rates of Ref. (l2| with those obtained from our 
calculations, leaving all other ingredients unchanged (We 
also change rates at iV=84 and 86, by amounts equal to 
the change in corresponding nuclei with iV=82.) By spec- 
ifying an appropriate temperature and density depen- 
dence on time, we mock up conditions in the "neutrino- 
driven wind" from Type II supernovae, the current best 
guess for the r-process site. 

The results appear in Fig. ^. As expected, the 



A « 130 peak shrinks noticeably. The A sa 195 peak 
broadens with the new half-lives because abundances 
around N—12Q are built up not just at the longest lived 
(most stable) nucleus produced, but at more neutron-rich 
A^=126 nuclei as well. As a result, more nuclei are popu- 
lated and the peak widens. By contrast, the N='bQ peak 
(not shown in Fig. |l^) does not change much; its shape 
depends largely on the half life of ^*Ni, which in our 
calculations is almost the same as in those of Ref. [|2|. 
We have already pointed out, however, that both life- 
times are probably too long. It is therefore reasonable 
to expect larger changes at low A than our simulation 
indicates. 

Besides altering the distribution of abundances, 
smaller half-lives can shorten the time required for the 
r-process to synthesize all the elements. The process 
proceeds only as fast as material can move through "bot- 
tlenecks", the especially long-lived isotopes at the three 
closed neutron shells. In many simulations the sum of the 
lifetimes of the bottleneck nuclei exceeds the expected 
duration of r-process conditions in neutrino-driven winds 
p6[ . To see what our lifetimes do, we run a series of r- 
process simulations at different temperatures, with both 
the half lives of Ref. [Of and with those presented here. 
At low temperatures, for which neutron photodissocia- 
tion rates are slower, the average nucleus is extremely 
neutron rich and so the change in beta half-lives does 
not have a large effect; the time required for the r-process 
drops by « 15%. At higher temperatures, nuclei that are 
slightly less neutron-rich are produced, and our half-lives 
have a more significant impact, resulting in an r-process 
time about 50% shorter. These high-temperature simula- 
tions have difficulty reproducing the observed abundance 
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FIG. 10. Predicted abundances in a simulation ol the 
r-process. The solid line corresponds to the rates of Ref. 
[nd, and the dotted line to the rates obtained here around 
N=82 and 126. All other nuclear and astrophysical parame- 
ters are the same for the two lines. The crosses are observed 
solar-system abundances. 



distribution, so we cannot at present take them very se- 
riously. It might be possible, however, to alter other 
conditions so that the correct distribution is restored. In 
any case, quantitative insight will have to await more 
comprehensive calculations of beta decay. 



V. SUMMARY AND CONCLUSIONS 



takes to complete the r-process will shrink, perhaps sub- 
stantially. But in order to explore these issues fully, one 
needs to know half-lives for all the waiting point nu- 
clei, not just those at the closed shells. Other quantities 
that affect abundance distributions, particularly neutron 
separation energies, also have to be better understood 



One virtue of our self-consistent framework is that 
many extensions and improvements can be made system- 
atically. First, a new interaction, based on the concept 
of the energy density functional, can be developed. This 
force should be able to reproduce global nuclear observ- 
ables, and at the same time yield a value oi g'g that is close 
to 1.8. Finding such a force will require abandoning some 
of the conventional relations between the Tig™'^ and Hg^'^ 
components of the energy density p4| ] . Then the parame- 
ters of Hg™" can be adjusted to global nuclear properties 
while those of Hg^'^ can be fit to spin-dependent proper- 
ties (go, ^0, moments of inertia, etc.). 

Following the development of a better interaction, one 
can both improve the calculations presented here and ex- 
tend them to open-shell nuclei that are not spherical. The 
formalism might also to be developed so that particle- 
number and isospin conservation are at least partly re- 
stored. Finally, all the nuclear ingredients in r-process 
network simulations, including neutron-separation ener- 
gies, should be based on the same effective Hamiltonian. 
Only then might our understanding of neutron rich nuclei 
be sufficient to predict details of the r-process abundance 
distribution. This paper is a first step towards that goal. 



This work contains the first consistent application of 
the coordinate-space HFB+QRPA formalism to beta de- 
cay in heavy neutron-rich nuclei. Our approach fully ac- 
counts for the coupling between the nuclear mean field, 
pairing, and the particle continuum. In addition our re- 
sults are based on one Hamiltonian; we use the same 
interaction, SkO', in the HFB and the QRPA calcula- 
tions. 

Beta-decay half-lives depend on Q- values, shell struc- 
ture, and the residual interaction. Much of the effect of 
the last ingredient is summarized by the Landau-Migdal 
parameter g'g. Apparently all Skyrme interactions com- 
monly used in nuclear structure studies, including SkO', 
have values of ^q that fall well below the experimental 
estimate (gQ^l.S). As a result, the centroid of the GT 
strength distribution is usually too low. That this de- 
fect is not fatal is due to pairing. Like-nucleon pairing 
produces diffuse Fermi surfaces, which in turn allow T=0 
pairing to pull strength down in energy. We can therefore 
compensate for a small (^q by weakening the strength of 
the T=0 pairing interaction slightly . But we definitely 
cannot omit it altogether; as Fig. ^ shows, calculations 
that do almost always overestimate lifetimes. 

As for the r-process itself, shorter beta-decay rates will 
have an impact, but just how much is not yet clear. 
Abundances certainly shift somewhat, and the time it 
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